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Abstract 

We examine the question of whether Bethe's ansatz reproduces all states in the periodic Heisen- 
berg XXZ and XXX spin chains. As was known to Bethe himself, there are states for which the 
Bethe momenta k n diverge: these are in fact the simplest examples of "string" solutions. The coef- 
ficients of the Bethe wavefunction, too, diverge. When there are only two down spins in the system 
(the case considered by Bethe), we can renormalize these coefficients to get a sensible (and correct) 
wavefunction. We show that this is not always possible when there are more than two down spins. 
The Bethe equations have several such divergent solutions, and some of these correspond to genuine 
eigenfunctions of the Hamiltonian, but several do not. Nor do they reproduce the correct energy 
eigenvalues. Moreover, we point out that the algebraic Bethe ansatz, an alternative way to construct 
the wavefunctions proposed by Faddeev, Takhtajan et al., leads to vanishing wavefunctions for all 
these solutions. Thus, the Bethe ansatz solution of the Heisenberg model must be regarded as either 
incomplete, or inaccurate. 

1 Introduction 

The XXZ spin-half Heisenberg Hamiltonian for N interacting spins, on a one-dimensional chain with 
periodic boundary conditions (n + N = n), is 

JV 

H = -2 J2 (S X nS: +1 + SlSl +1 + AS* n S* +1 ) (1) 

n=l 

where A is the anisotropy parameter, for the XY model, 1 for the isotropic ferromagnet. A = — 1 can 
be mapped to the isotropic antiferromagnet. We have put the overall factor of 2 for convenience. In 
terms of the Pauli spin matrices, 

x • A 



H = - J2[ a n °n+l + °n °n+l + ( 2 ) 



2 

where a ± = (a x ± ia y )/2. This was among the first many-body quantum problems to be solved exactly. 
Bethe [|J solved the isotropic model (the "XXX model", A = 1) in 1931, and the XXZ Hamiltonian can 
be solved similarly with little additional effort. The form of wavefunction which he used has become 
known as the Bethe Ansatz and has been applied to a wide variety of other problems. 

The Bethe ansatz, however, is a guess, and while one can generally verify whether or not it works, 
there is no easy way to tell whether it is "complete" , that is, whether all the energy eigenstates are of 
this form. 
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There are actually two versions of the Bethe ansatz — the original method of Bethe, and the "algebraic 
Bethe ansatz" developed by Takhtajan, Faddeev and others 0, || based on the work of Baxter (J] on 
vertex models. Though they start from very different ideas, the two approaches lead to the same equations 
and (apart from normalization) the same wavefunctions. Some concrete statements on completeness have 
been made in the case of the XXX model, namely that the algebraic Bethe ansatz reproduces only the 
states of "highest weight" with respect to the underlying SU(2) symmetry (which has been proved) and 
that it reproduces all such states (which has not been shown strictly). By "highest weight" is meant 
a state with the maximum allowed value of S z (the total z component of the spin) for given total spin 
S; such a state is annihilated by the total spin raising operator S + . The proof that all Bethe states 
are of highest weight ||, however, depends on a change of variables from Bcthe's k n to new variables 
A„ which arise more naturally in the algebraic method (besides being convenient in other ways), and 
the subsequent restriction that these variables must be finite- valued and distinct. The original ansatz of 
Bethe, in terms of the fc's, reproduces many non-highest-weight states without trouble. Moreover, this 
doesn't answer the question of completeness in the anisotropic (XXZ) case where there is no underlying 
SU(2) symmetry. The arguments for completeness even in the XXX chain are not rigorous: Essler et 
al. p], for instance, show that the usual state-counting argument [|| |3| is not quite correct, but they 
view this as a reorganization of states and suggest that the ansatz for the XXX chain is SU(2) complete 
anyway. 

Here, after a quick review of the Bethe ansatz, coordinate and algebraic, mainly to fix the notation, 
we look into these questions. We find that a number of states correspond to Bethe fe's which diverge 
to ±ioo, and these include states which are of highest weight at the isotropic point. The variables A„ 
which one usually uses in the algebraic Bethe ansatz are well-behaved, but the coefficients of the Bethe 
wavefunction are singular. The algebraic method offers an alternative way to produce a wavefunction, 
but this construction produces a vanishing wavefunction. 

To handle this situation, one can consider only the ratios of the Bethe coefficients, rather than the 
coefficients themselves (which is reasonable since the wavefunctions are not normalized) . These turn out 
to be finite, but we then find that some of the states thus produced are not eigenstates of the Hamiltonian 
at all. In fact, the energies predicted by the Bethe equations also turn out to be wrong. 

In short, the equations of Bethe, written in terms of the variables A„, have a number of solutions with 
finite, well-behaved A n which (depending on one's viewpoint) either do not reproduce the corresponding 
eigenstates, or produce too many states including some which are not eigenstates. 

Finally, we also point out that one can formally get rid of these singularities by introducing an 
Aharonov-Bohm flux through the ring, that is, associating different phases with the forward and backward 
spin hopping amplitudes. But for certain special values of the phase (namely, the N roots of unity — 
corresponding to a total flux of a multiple of 2ir through the ring) the singularities reappear, and for values 
of the flux close to zero or to any of these points the Bethe equations become numerically ill-behaved and 
any rootfinder would have severe problems converging to a solution. 



2 The coordinate Bethe ansatz 

Since the Hamiltonian (Q) commutes both with the total spin and with a z = Y] n a z t , we can work with 
an Ising basis with fixed a z . Suppose the number of down-spins is I (< N/2: the other states can be 
reached by symmetry), and for a given Ising state, suppose the positions of the down-spins are xx, x%, 
. . . , xi, and call this basis state \x1X2 . . . x{). A general wavefunction can then be written in this basis as 

I*) = X! ip{xi,x 2 ,...,xi)\xx,X2,...,xi). (3) 

X 1 <X 2 ...<Xi 

Bethe suggested the following form for the expansion coefficients, 

ij>(xi,X2, ...,xi) = ^2 ^(-P)exp i^2k Pn x n J (4) 
p \ n=l / 

where P is a permutation of the integers 1, 2, . . . , I, the sum is over all permutations with amplitudes 
A(P), and the fe's are some as yet undetermined quantities. 
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One can show that such a wavefunction solves the Hamiltonian (0) if the amplitude for exchange of 
two neighbouring particles is 

_ A(...,m,n, ...) _ e i( fc m+ fc n) - 2Ae lfc ™ + 1 
mn = A(...,n,m,...) ~ ~ e *( fe ™+^) - 2Ae* fc ™ + 1 ' ( "* 

In addition, periodic boundary conditions imply that 

e — = n *™ = (-i)'- 1 n l^i'lk+i - (6) 

Equations @ can be solved for fc n , and each solution {fc„} gives us an eigenstate of (g), with associated 
energy 

I 

£ = -7VA/2 + 2ZA-2^cos(fc„). (7) 



n=l 

For notational simplicity, we will often use 

so that the Bethe equations become 

n / - \ /. i i i z m z n 2Az m -(- 1 



(8) 



*»" = (-I)'" 1 II _a It (9) 
3 The Algebraic Bethe Ansatz 

The same equations can be derived rather differently, by diagonalizing the transfer matrix of the "six 
vertex model" , in which the Hcisenberg XXZ Hamiltonian is embedded. Here we will concern ourselves 
only with the XXX model, where the algebra is much simpler. We will not describe this method (see ||, 
for instance) but only mention the change of variables involved, and how the eigenvectors are constructed. 
In the case of the XXX model, we define new variables A„ as follows: 

An + »/2 nnA 
Z - = X—-f2 (10) 

A„ = l - Z -^± = icot(fc„/2). (11) 

These variables arise very naturally in the algebraic approach. In terms of these, the Bethe equations (|J) 
become 

^±m N = n (12) 

\ m — i 2 J V- \ m — \ n —i 

n^ra 

and the energy of the corresponding state is 

N 



E = - 



— (13) 
2 ^1 + 4A„ 2 



The commuting family of transfer matrices for the six vertex model, T(A), embed the Heisenberg 

Hamiltonian, and the eigenvectors of the former are the eigenvectors of the latter. T(A) is the trace of 
the monodromy matrix 

r(A) =i w (A)Ljv-i(A)...ii(A) (14) 
where L n (X) is the "local L operator" 

Ln(X)=( A + (i/ + 2) < , z V (15) 

V l < A-(i/2K J V 1 
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Writing the monodromy matrix as 



it can be shown that the vector 

|{A„» = {j[ B ( x n)j im---> (17) 

is an eigenvector of the transfer matrix T, and therefore of H, if the {A„} satisfy the equations (|T^). 

4 Completeness 

Since the Hamiltonian (||) is translationally invariant, a useful classification of the energy eigenstates is 
in terms of eigenstates of the translation operator T which shifts all spins by one lattice site to the left. 
All eigenstates of H can be chosen to be eigenstates of T (non-degenerate eigenstates are necessarily so 
anyway) with eigenvalues 4> such that cf> N — 1. The Bethe ansatz does this automatically: for given {z n }, 
(/) = l\z n . 

There is a limit where the Bethe method does reproduce all states without trouble: that of A = 0, 
or the XY model, which is equivalent to a free Fermi system. The exchange amplitudes A mn in (^|) each 
become —1, and the fc's are 2ni/N times an integer (for odd I) or a half-integer (for even /), all distinct 
(otherwise the wavefunction vanishes). We then have exactly as many states as we need. Some of these 
states evolve smoothly as one turns on A, but some don't, and it's interesting to see why not. 

It turns out that the A = states are sometimes degenerate, but the degeneracy is lifted for any finite 
A. The states, if chosen correctly, do evolve continuously as A is changed from zero, but the A = Bethe 
ansatz gives superpositions of these which are no longer eigenstates for any finite A. We see concrete 
examples of this below. This also suggests a way around the problem: introduce an Aharonov-Bohm 
flux, which does not change the eigenstates at A = but lifts the degeneracies in energy. We discuss this 
later. 

4.1 Four sites 

The generic Bethe Ansatz wavefunction for two down-spins is 

M = E (V^2 X2 + ^Z2 X1 Z1 X >) \ Xl ,X 2 ) (18) 

with x\ < X2] and we also have 



3=1, X2 



zi 4 = A, (19) 

z 2 4 = 1/A, (20) 

A = _*i*2-2Azi + l_ (21) 
ziz 2 — 2Az 2 + 1 

We will number sites starting from 0, to ease the algebra. 

Let us follow the convention of denoting a down-spin by 1 and an up-spin by 0. In the case of four 
sites and two down-spins, we can write down the eigenstates of T, with corresponding eigenvalues </>, 
immediately: 





= |1100) + 


|0110) - 


h |0011) - 


f |1001) 




= |1010) + 


|0101) 








= |1100) + 


i|0110) 


- |0011) 


— ijlOOl) 


liM 


= |1100) — 


|0110) - 


- |0011) 


- |1001) 




= |1010) — 


|0101) 








= |1100) — 


i|0110) 


- |0011) 


+ i|1001) 
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Now, for (p — one has to further diagonalizc H in the subspace of and l^j); but it turns out 
that all the rest are already eigenstates of H . Moreover, ip4 and ips at A = are degenerate, and the 
Bethe ansatz gives us some superposition of these; but the degeneracy is lifted for arbitrarily small finite 
A. This is one reason the completeness question is nontrivial: it is not obvious that new well-behaved 
solutions of the Bethe equations will emerge for finite A, and if they do, they must be ill-behaved as 
A — > since we already have all the solutions there. 

In fact, ^5) is reproduced with the choice (2:1,2:2) = ±1- This solution works for all A except 0. We 
now demonstrate that there is no well-behaved solution for {tpi}- 

From the generic Bethe wavefunction in this case (fl8|), we see that Z1Z2 = —1- On the other hand, 
the coefficient of 1 1010) should be zero, so we have Z2 2 + (l/A)zi 2 = 0. Combining these conditions gives 
us zi 4 = ~A which contradicts ( |l9| ) unless z\ = 0, which implies A = 0, z-i = 00. 

All the other states evolve smoothly as A is turned up from zero; however, singularities occur at 
particular values for A. At A = 1, the upper of the two <fi = states corresponds to z% = Z2 = 1 so 
the Bethe wavefunction is singular (or zero). The same thing happens to the <fi = ±i states at A = \/2. 
Generally this "collision" of z's or fc's occurs at the point where the fc's go complex. One can make sense 
of it, however, by considering the limit of the wavefunction as A approaches the critical value. 

In terms of the algebraic Bethe ansatz (we confine ourselves here to the XXX point), things are a 
little different: if we require that the A's, rather than the /c's or the z's, be finite, the value z = 1 is 
forbidden. Then the states \ip3), \ip5) and \tpe) (which are not highest weight states) are not reproduced. 
It may appear that the situation is saved for the highest-weight state \ip4): the values and 00 for z 
translate to the nicer values ±i/2 for A. The energy ( |l3| ) appears to be divergent, but if we take the A's 
to be (—i/2 + e, i/2 + e) and let e vanish, we do get the correct value. Unfortunately, these A values don't 
produce an eigenstate: B(i/2)B(— i/2)\ TTTT) = 0, which is easy to verify. 

So there exists at least one state, \ip4), which is not reproduced by the algebraic Bethe ansatz at the 
isotropic point. In fact a similar singularity exists with the state 

£(-i)V-a- +1 |TTT--->, 

for all even N. This state was actually known to Bethe H), who remarks that the fc's must diverge to 
accommodate it. 

4.2 Six sites 

Here things get a little more complicated and we have many more states to contend with. But classification 
of states in terms of (f>, the eigenvalue of the translation operator T, can still be done. In case of a real 
translation eigenvalue <f> (±1) we can further consider the parity operator, which reverses the order of the 
spins. (If 4> is not real, this doesn't commute with T, but one can generalize the definition to take care 
of that.) So consider the state 

= 1110100) + 1011010) + 1001101) + . . . 

-|101100) - 1010110} - ... (22) 

which is the only state with 0=1 and parity —1, and is therefore an eigenstate of H for arbitrary A 
(also, as one can verify, a highest weight state). Consider also the state 

|V/) = 1110100) - 1011010) + 1001101) + . . . 

+ 1101100) - 1010110) + 1001011) - . . . (23) 

which is also an eigenstate of H (the only one with <p = — 1 and parity 1), and a lowering of a state 
analogous to ^4) in the previous section. Neither of these states can be reproduced unless either two of 
the z's are equal, or one of them vanishes (and, therefore, another diverges). Since the proof is not as 
compact as in the four site case, we leave it for appendix These states, too, are degenerate with other 
states at the XY point, and some superpositions of these are reproduced by the Bethe ansatz. 

In appendix |^, we present the values of A which produce these states at the isotropic point if we take 
appropriate limits for the divergent Bethe coefficients. As predicted, the corresponding values of k are 
divergent. 
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4.3 Adding a flux 



The XY degeneracies can be lifted by passing an Aharonov-Bohm flux through the ring. This modifies 
the forward and backward hopping amplitudes by phase factors: 



N 



H = }_,[ "^n °n+l + aG n Vn+1 + ^ «+l ) > ( 24 ) 
n=l ^ ' 

where a is a complex number of unit modulus. For a = 1 we recover the original Hamiltonian. For most 
other values of a, the ground state degeneracies are lifted and the states evolve smoothly when one turns 
on A. But the singularities reappear when one removes the flux. Moreover, certain other values of a also 
prove to be singular. This is an example of the fact, noted by Byers and Yang JtJ, that the energy levels 
of a quantum system on a ring are periodic in the flux through the ring. To see it in this case, we carry 
out a Bethe ansatz solution and get the analogue of (||): 

Zm = ( _1 J II ~ — i — ' ( 25 ) 



a*z m z n - 2Az n + i 



which, if we define w n = az n , becomes 

JV„ N 



Wn 



a N = (-l) i_1 TT WmWn ~ 2AnJm + 1 (26) 

*J- W m W n - 2Aw n + 1 



These are identical to the Bethe equations if a N = 1, so for these values of a the singularities will persist. 

The exact solution of the <j> = — 1 sector for four sites with flux is given in appendix [b|, and its 
behaviour as the flux is made to vanish can be seen explicitly there. 



4.4 Longer chains 

The common thing to the above examples is the divergence of one of the z's, or of the k's to ±ioo. In fact 
for any system size and any number of overturned spins, if one of the z's diverges at A = 1, the algebraic 
Bethe ansatz fails to reproduce the corresponding state. To see this in the XXX case, note that if one z 
diverges, another must vanish; the A's corresponding to these are ±i/2. Consider the general algebraic 
Bethe state: 

\{K}) = l[B(x n )\m ■••> 

Choosing the argument of the last B operator to be i/2, we can easily convince ourselves (from the 
definition of B in terms of the L matrices in sec. [|) that the action of this on the ferromagnetic state is 
merely to flip the first spin. If we choose the next B operator to have the argument —i/2 that state is 
killed entirely. The values of the other A variables then don't matter. 

It is also quite clear that such solutions will in general exist for all lattice lengths. To see this, plug 
the values A; = i/2, A;_i = —i/2 into eq. (|l2]). These equations are immediately satisfied for rn — I and 
m = I — 1. For the rest, they now become 

(A m + A m ~ 3i/2 _ l r A m 

Am - i/2 J A™ + 3i/2 ~ W 

and any solution of these for Ai, . . . , X1-2 combined with the above values for A;_i, A; is a solution of 
©• 

This, by the way, is the simplest sort of "string" solution — such solutions, of groups of complex A n 
with a common real part and the imaginary parts separated by i, arranged symmetrically about the real 
axis, are generic in the thermodynamic limit. This particular string exists for all system sizes. 

Therefore for any lattice length N and spin sector I there exist several solutions for the equations ( |l2| ) 
for which the algebraic Bethe ansatz gives a vanishing wavefunction and the coordinate Bethe ansatz a 
singular one. How to handle these singularities, and whether they can give us meaningful wavefunctions, 
are discussed later. 



A„ + i 



— A,, — i 



(27) 
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Though we have discussed only the XXX model above, the same is true of the XXZ model. The 
variables A„ are defined a little differently here (the earlier ones are related to these by a scale factor) 
and the Bethe equations are H 

sin(A m +r)) \ N _ -r-r sin(A m - A» + 2rj) 
sin(A m - rf) ) ±L s in(A m - A n - 2?;) ' 

Here A = cos2?7. These have the singular solution A/_i, A; = ±77 with I — 2 equations for the remaining 
A m . The L matrix in this case is 

j m _ ( W4(A) + w 3 (X)cr n z (sm2ri)a n " \ . . 

Ln[A > \ (sm2n)a n + w 4 (\) - w 3 (\)(Jn z J { ' 

with (1V4 + Ws)(X) = sin(A + 77) and (1V4 — W3XA) = sin(A — rj). The B operator is defined as before, and 
one can convince oneself that this choice of A„ annihilates the reference state, just as in the XXX case. 
Unfortunately our proofs regarding the specific singular states are not readily generalized to longer 



chains, so other than the one example at the end of sec. 4.f, and the states suggested at the end of 
appendix we have no rigorous examples to give. But the above examples involve states with 5- 
independent coefficients, which are degenerate with other states at S = 0, and such states certainly exist 
for larger lattices as well, so it is entirely conceivable that they too could pose problems for the Bethe 
ansatz. 

We can be a little more specific. In the XY limit, the fc's are the "bare" momenta 2TrI n /N, where I n 
are distinct integers (for odd I) or half-integers (for even I), Consider the sector <f> = 1 (cf) being, as above, 
the eigenvalue of T); this is the same as saying ^2 k n = 0. This can happen in two ways: either the k n are 
symmetrically distributed around 0, or they are not. If they are not, the solution corresponding to — k n 
is distinct from this solution but degenerate in energy with it. However, these two solutions, generally, 
do not evolve into well-behaved solutions as one turns on A: the solutions at finite A correspond to 
superpositions of these, which are no longer degenerate. The same problem occurs in the kn = tt 
sector. 

Numerically, we investigated several such states. The idea was to solve equations (||) using a Newton- 
Raphson rootfinder (capable of finding complex roots: a simple modification of the standard method). 
For this method to work, one needs a good starting guess. So a possible method is to keep a finite flux $, 
start the rootfinder from the XY point where the solution is known (and not affected, except in energy, 
by the flux), evolve A forward to some desired value and then evolve the flux down to zero. In every 
case, the last step failed to converge, yielding divergent z's or k's or else hitting some spurious solution 
which has two k's equal and therefore yields a vanishing wavefunction. 

It seems therefore that such singularities are very common, at least in the (j> = I and cj) — — 1 sectors, 
possibly in the other sectors too. 



5 Handling the singularities 

Thus, for several eigenstates of the Heisenberg model the Bethe ansatz becomes ill-behaved and the 
coefficients of the wavefunction, diverge (in the coordinate Bethe ansatz) or vanish (in the algebraic 
method). For the coordinate method, one can handle the divergence by considering only the ratios of 
the coefficients, rather than the coefficients themselves. For the algebraic method, since the problem 
values of A arc ±i/2, one can let them approach these values and consider the limiting ratios of the 
coefficients of the wavefunction. While this should work in principle, it would be rather messy and we 
don't pursue it further: the algebraic BA is valuable more as a formal elegant unifying tool than as a 
practical calculational aid. 

For the four site, t wo spin down case, if we take Ai = i/2, X2 = —i/2, this works fine and we reproduce 



the state ^4) of sec. This is demonstrated in appendix |B|. However, consider the six site case with 
three down spins. If we choose A2 and A3 to be H/2, we are left with one Bethe equation for Ai which 
turns out to be a fifth degree equation with six solutions if one includes Ai = 00 (or z\ = 1). That, and 
Ai = 0, do lead to valid eigenstates which are precisely the problematic eigenstates discussed in sec. |4.2|. 
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But the other four solutions for Ai, though perfectly valid solutions of the Bethe equations (|12]), do not 
lead to correct wavefunctions (nor, indeed, to correct energies). This is further discussed in appendix 
For any even lattice length with three down spins, the Bethe equations have N such singular solutions 
of which two yield genuine eigenstates and the rest are spurious. This presumably happens for all spin 
sectors, but we don't have good estimates on what fraction of the total number of states is so affected. 

Another possible way to handle these singularities (when solving the equations numerically) is to 
introduce a finite flux, and slowly turn it down to zero. It turns out that the equations are so ill-behaved 
that this does not help. A rootfmder such as the Newton-Raphson method, which is based on the 
Taylor series and requires smoothly varying functions, behaves very badly in the neighbourhood of such 
singularities. 

Other kinds of singular Bethe states are possible — for instance, at certain special values of A, two fc 
values may collide. If the corresponding Bethe state has a well-defined limit as A — > Ao, where Ao is the 
value where the fc's coincide, then there is no problem. This happens for instance with the <f> = ±i and 



one of the cf> — 1 states in sec. iA t at S = y2 and 5 = 1 respectively. And then again, the fc's may be 
well-behaved but the A's singular, as in the state 1 1010) — |0101). Here, if one just works with the fc's, no 
problems arise. In particular, the Bethe coefficients are well-behaved. 



6 Conclusion 

The Bethe ansatz is a particular guess for a wavefunction and substituting it in the Hamiltonian, and 
applying periodic boundary conditions, leads to certain conditions on the parameters of the wavefunction, 
the fc's or the A's as one may prefer. We find that though Bethe's equations for the A's by themselves have 
a large number of solutions, some of these solutions lead to singular, ill-defined wavefunctions. These 
problems arise from divergences in Bethe's original quasimomenta fc„, though the A n to which they can 
be transformed are well behaved. 

We give explicit examples of this for four and six site lattices, but show that such solutions of the 
Bethe equations should occur for all lattice sizes. The question then is whether to allow such solutions 
or not. If we do not allow them, the Bethe ansatz is incomplete. If we do allow them, we get vanishing 
wavefunctions with the algebraic Bethe ansatz, and divergent ones with the coordinate method. If we 
attempt to make sense of these by considering only ratios of the Bethe coefficients, we get well-defined 
wavefunctions but not all of these are eigenstates of the Hamiltonian. Indeed, even without calculating 
the wavefunctions, we can check that the energies are wrong. 

Earlier arguments for completeness have been based on counting of solutions to the Bethe equations 
(|l2|), using the string hypothesis. But as we have seen, not every solution to the Bethe equations 
produces an eigenstate of the Hamiltonian. Therefore, even if the Bethe equations themselves have as 
many solutions as we require, it does not imply completeness of the method. 

In the case of two down-spins, however, the Bethe ansatz is complete if one normalizes the divergent 
wavefunctions appropriately. This was known to Bethe as was the fact that the fc's and the amplitudes 
for the wavefunction diverge for one state in this case. There are no incorrect solutions in this case. 

One question to be answered is how many incorrect solutions do we obtain in the thermodynamic 
limit. In the case of three down spins, we get N — 2 wrong answers, which is a vanishing fraction of the 
total number of states in the thermodynamic limit. It will probably remain a vanishing fraction if the 
number of down spins is finite, but if a finite fraction of the spins are down — if S z = 0, for instance — we 
haven't formed an estimate. 

One way to get around this, in principle, is to add a small flux to the system; but this helps only in a 
formal sense. The equations are very ill behaved near these singularities and any rootfmder would have 
a great deal of trouble converging to a solution. 

Finally, some remarks on why all this is important. One reason is a matter of principle it is good 
to know to what extent a given method solves a problem. Another is that the question could be relevant 
to other problems too — for instance, the Hubbard model, which becomes the Heisenberg antiferromagnet 
in one limit, may also exhibit similar singularities in its Bethe ansatz solution. This incompleteness will 
probably not matter much in calculations of thermodynamic properties, since the singular states may 
then be a vanishing fraction of the total number of states. However, when considering small systems, 
or systems with very few overturned spins such as low-lying excitations above the ferromagnetic state, 
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perhaps one should keep the completeness question in mind. And finally, the above considerations are 
very important indeed if one wants to use the Bethe ansatz for numerical calculations on a finite sized 
lattice. This should not be a common difficulty since other methods, such as exact diagonalization, can 
easily be used in one dimension. 
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A Singular states in the six site lattice 

For six sites and three down-spins, we have three z variables and three two-particle scattering amplitudes 
A mn (since A nm = l/A nm . Let A X23 = 1; then we have 

A _ 1 A _ 1 A _ 1 

^132 — —. , -A213 — —. , ^231 — 



^23' ^12' ^12^13' 

A 312 = X~~ 4~ ' A 32i = A l A - (30) 

^23^13 ^-23^13^12 

Also, from the Bethe ansatz equations we have 

Zl 6 = A 12 A 13 , z 2 6 = ^-, zs= 1 r^r- ( 31 ) 

A\2 ^13^23 

Let us define A 23 = A and rewrite A V2l A 13 in terms of this and the z's. Using all these in the formula 
(Q) for the Bethe wavefunction, we get for tp(xi, x 2 ,x 3 ) the formula 

il>(xi,X2,x 3 ) = z 1 ^z 2 X2 z 3 X3 + jz 1 ^z 3 ^z 2 X3 + jz 2 ^ +e z 1 ^z 3 X3 

+z 2 x 'z 3 X2 Zl X3 - 6 + z 3 Xl+6 z 1 X2 z 2 X3 + -z^z^z^ (32) 



Now consider the state ( |22| ) in sec. 4.2. To reproduce this state, we require that 

^(0,1,2) = 0, (33) 
^(0,2,4) = 0. (34) 

Each of these, together with Z\z 2 z 3 — 1, can be used to solve for A. Using we get 

z 2 2 

A = -A, (35) 

zr 



and using (33), 



zi 3 z 2 + Zl 3 Z 2 4 + z 2 



A = — -T^l ~ (36) 

Zl 3 Z 3 + Zl 3 Z 3 4 + Z 3 

Unless z 2 or z 3 is zero, we can equate these, and after cancelling a common z 2 /z 3 factor, rearranging and 
using z\z 2 z 3 = 1, we get 

zi 3 (z 2 - z 3 ) - Zl 2 (z 2 2 - z 3 2 ) + {z 2 -z 3 )=0 (37) 
and if we multiply the last term by z\Z 2 z 3 , this can be factorized as 

zi(z 2 - z 3 )(z 3 - z\){z\ - z 2 ) = 0. (38) 

Thus, if the variables z\, z 2 , z 3 are all distinct, one of them must vanish (and consequently another 
must diverge). 

A precisely analogous proof goes through for the state (p3|) using the condition z\z 2 z 3 = — 1. 
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B Some exact solutions in the four site case 



Here, as in section [|, we work with states labelled by 4>, the eigenvalue of the translation operator T 
which shifts all spins by one lattice site to the left. 

B.l Sector <fi = — 1 with flux a 

We have two solutions in this sector, given by 



-7 ± y 7 2 + 4 

(zi, z 2 ) = \ (39) 



where 7 has the two possible values 



7 - ±±^3. (40) 

A - (41) 
a — a* 

Note that as a — > 1, A diverges, and consequently one solution for 21,% diverges, while the other 
approaches ±1. 

B.2 Sector <f) = — 1 with no flux 

The state 1 1010) — |0101) is readily reproduced with z = ±1. To reproduce the other state, we need 
z\ — » 0, Z2 — > 00. Now the Bethe coefficient of the state 1 1010) is Z2 2 + z\ 2 /A, which on using Z1Z2 = — 1 
and zf = A becomes 2z2 2 - The Bethe coefficient of 1 1100), similarly, is z^ — Z2 3 - Both of these diverge 
as Z2 — > 00, but if we take their ratio, this vanishes. (It also vanishes, as it should, if Z2 — > 0.) So for 
this state, we can conclude that the 1 1010) basis state is of zero weight compared to 1 1 100) , and we do 
reproduce the correct eigenstate. 

B.3 Sector 0=1 with no flux 

There are two solutions in this sector, given by 



= 7±V f^ (42) 



for each of the two values of 7 



A±VA^+8 

7= g ■ ( 43 ) 

As A — * 1, one of the 7 solutions tends to 2, for which the z values collide at the value 1. At this point, 

Z1Z2 — 2Azi + 1 



.4 



ziz 2 - 2Az 2 + 1 



also becomes ill-defined. If it weren't so, the wavefunction would vanish. Now, as before, we can consider 
the coefficients of 1 1010) and 1 1100) and eliminate the ill-behaved A from them, before taking the A — > 1 
limit. The former, as before, can be written solely in terms of Z2 as 2z2 2 and the latter as Z2 + Z2 3 - Both 
of these have well defined limits as Z2 — > 1, and the correct wavefunction is reproduced. 

C Six sites: exact solutions for the singular states 

Here again we confine ourselves to the XXX model. Let us look for a solution of the equations for A m , 
(|i"2|), for six sites and three down spins, with Ai = — i/2 and A 3 = i/2. Then we only have to find A 2 , 
which we will call A. Equations (|l^) are already satisfied for Ai and A 3 , while for A 2 = A they become 



A + i/2\ 6 _ f\ + 3i/2\ ( X + i/2 



A -i/2/ V A -i/2 / \\-3i/2 



(44) 
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which have solutions A = oo (which we will not disallow), A = and the four roots of 



16A 4 + 40A 2 -7 = 0. (45) 

For now, we will not specify which root we are considering; for all of these, the Bethe equations (|l^) are 
satisfied. The z values corresponding to these are Z\ = 0, 2:3 = 00 with various finite values of z-i. +1 for 
A = 00, —1 for A = 0, and other values for the other four solutions of A. 

Our Bethe wavefunction ( p2| ) is then singular because of the singular values of the z's. Let us, however, 
assume 

{z 1 ,z 2 ,z 3 ) = (l/(3,z,ct>P/z) (46) 

and take the limit of this as (3 — > 00. Here <f> is, as before, the eigenvalue of the translation operator: 
cf) 6 = 1. In that case, to leading order in (3, the coefficients of various Ising states in the Bethe wavefunction 
are 

1 



|11 1000) 




1110100) 




1101100) 


z z 


(101010) 


,4,1 + 



with coefficients of other Ising states being the above multiplied by appropriate powers of Z1Z2Z3. As 
(3 — > 00, the coefficient of 1 101010) becomes insignificant. First consider the case z-i = z = ±1, or 
A = 0, 00. If we choose <fi = ±1 as appropriate, the coefficient of 1 111000) can also be made to vanish (to 
order f3 5 ) and we recover precisely the two states discussed in appendix 

But now consider the other four possible values of z. In all eigenstates, the coefficients of 1 101 100) 
and 1 110100) differ by a factor of <f> (or a power of it): (f>, remember, is a sixth root of unity. Here they 
differ by a factor <fi 4 /z, which is not a power of (f> since the four values of z other that ±1 are not sixth 
roots of unity. So for these values of z, the states (which are in any case singular) are not eigenstates of 
the Hamiltonian even in a limiting sense. 

In fact we have checked the parentage of all the exact eigenstates in the six site case, and apart from 
the two discussed above, all are reproduced with finite z n , with at worst a collision of two values. So 
there is no room for the remaining four solutions. Moreover, if one uses the values of A calculated in this 
appendix to calculate the energy by equation (|l^), using a limiting method as in sec. 4.1 to handle the 



divergence for Ai, A3 = ±2, one gets the correct answers for A2 = 0, 00 but wrong answers for the other 
four solutions for A2 — as one can verify with exact diagonalization. (For this calculation, all we need to 
assume is that 

4 4 

= 2, (47) 



l + 4(i/2) 2 1 + 4H/2) 2 



which as a limit is justified in sec. 4.1, and which produces the correct energy when the wavefunction 



is correct). So there exist solutions of the Bethe equations (12) which do not correspond to actual 
eigenstates of the problem. 

In the case of three down spins, the above generalizes quite readily to any even number of sites: for N 
sites, we get N — 2 spurious solutions, and two solutions for A which correspond to genuine eigenstates. 



One of these is a lowering of the two down spin state given at the end of sec. 4.1; both remain eigenstates 
for all A. Thus, for eight sites, we get for (j> = 1 

l^i) = 111010000) - 1 10110000) - 111001000) + 110011000) + translations with = 1 (48) 

and for <f> = — 1 

^2) = 111010000) + 110110000) + 111001000) - 110011000) + translations with cj) = -1. (49) 
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